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We analyze two-color spatially localized modes formed by parametrically coupled fundamental 
and second-harmonic fields excited at quadratic (or x' 2 ') nonlinear interfaces embedded into a 
linear layered structure — a quasi- one- dimensional quadratic nonlinear photonic crystal. For a 
periodic lattice of nonlinear interfaces, we derive an effective discrete model for the amplitudes of 
the fundamental and second-harmonic waves at the interfaces (the so-called discrete \^ equations) , 
and find, numerically and analytically, the spatially localized solutions — discrete gap solitons. For 
a single nonlinear interface in a linear superlattice, we study the properties of two-color localized 
modes, and describe both similarities and differences with quadratic solitons in homogeneous media. 



I. INTRODUCTION 



The physics and applications of photonic band gap ma- 
terials (or photonic crystals) have been an active topic of 
research for more than a decade jjj. The next step in 
application of photonic crystals is to create tunable band 
gap materials where the gap could be controlled by an 
external parameter. One of the recent suggestions j^] 
is based on the inverse opal structure — a microscopic 
lattice of spheres of air surrounded by silicon — with a 
layer of liquid crystal material that can make the trans- 
mission properties programmable by applying an electric 
field. Continuous tuning can also be realized by varying 
the liquid crystal temperature || . 

Another important idea towards creating dynamically 
tunable band gap materials for switches and transistors 
operating entirely with light is to employ their nonlin- 
ear properties, thus creating nonlinear photonic crystals. 
The concept of nonlinear photonic crystals, defined as 
having a spatially periodic nonlinearity, was introduced 
by Scalora et. al. J4j in a numerical study of ultra-fast 
optical switching and limiting in cubic nonlinear Kerr 
materials. Bistability and localized modes in photonic 
superlattices with embedded layers possessing nonreso- 
nant cubic (or Kerr-type) nonlinearities have also been 
discussed in the literature ||. Recent advances in the 
so-called cascaded nonlinearities Q demonstrate an ef- 
fective way to lower the switching power by employing 
parametric interaction and frequency conversion in non- 
centrosymmetric quadratic nonlinear optical materials. 
Parametric interactions are also known to support soli- 
tary waves, spatial quadratic solitons Q, that exist in 
homogeneous media where spatial localization is induced 
by two-wave parametric mixing processes between the 
fundamental wave and its second harmonic. Quadratic 
nonlinear photonic crystal was introduced as a concept 
by Berger || in a study of multi-wavelength frequency 
conversion. 

Taking into account the similarities between the lo- 



calized defect modes in linear inhomogeneous media and 
nonlinear localized modes in homogeneous media [0], we 
wonder if parametric interactions can support localized 
modes in a variety of photonic band gap structures with 
quadratic (or x^) nonlinearities. The first step in this 
direction has been recently presented in Ref. ML where 
we have analyzed second-harmonic generation (SHG) at 
a thin, effectively quadratic nonlinear layer separating 
two (generally different) homogeneous linear media, and 
predicted multistability of SHG for both plane waves and 
localized modes, also describing two-color localized pho- 
tonic modes that can be excited at the interface. 

The main purpose of this paper is twofold. First, we 
generalize the results of Ref. |1(J to the case of a thin 
quadratic nonlinear layer embedded into an arbitrarily 
stratified periodic linear medium. In particular, we con- 
sider a nonlinear defect layer with second-order nonlinear 
response in a perfectly periodic dielectric structure — a 
one-dimensional analog of a photonic crystal with a non- 
linear impurity. Secondly, we develop a general formalism 
for analyzing localized modes in multilayer structures - 
nonlinear superlattices — and describe two-color localized 
gap modes supported by a periodic lattice of thin layers 
with quadratic nonlinearity. 

The structure of the paper is the following. In Sec. || 
we present our model that is described by a system of two 
coupled nonlinear equations for the envelopes of the fun- 
damental and second-harmonic waves. In the stationary 
case, the solution can be presented as a superposition 
of the forward and backward traveling waves, and this 
allows to derive an effective system of discrete coupled- 
mode equations for the wave amplitudes at the layers, the 
so-called discrete x*^ equations (Sec. Ill A). Solutions of 
these equations for localized nonlinear modes are briefly 
discussed in Sec. [II B. At last, in Secjfv] we consider a 
generalization of the results of Ref. [flCfl to the case of a 
single quadratic nonlinear layer embedded into a periodic 
linear medium, and describe a number of new features of 
the localized modes that appear due to a bandgap struc- 
ture. 
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II. MODEL 



First, we discuss the physical motivation for our model. 
Let us consider an interface between two semi-infinite 
bulk optical media with inversion symmetry. The inter- 
face layer breaks the symmetry and therefore it should 
possess a nonvanishing surface quadratic response Jll]] 
that can be enhanced by a proper coating creating a 
nonlinear layer with quadratic nonlinearity |T^ |. Such 
a layer corresponds to an effective nonlinear defect that 
can support photonic modes localized at the interface. 
There exists a strong experimental evidence of SHG in 
such type of localized photonic modes. For example, re- 
cent experimental results jl3) reported SHG in periodic 
photonic band-gap structures with embedded nonlinear 
defect layers. An enhancement of the parametric interac- 
tion in the vicinity of the defects was observed, suggesting 
that SHG occurs in localized modes, being suppressed for 
propagating modes. The surface nature of SHG was con- 
firmed by comparison of experimental data and results 
of direct numerical simulations jL2j . 

To introduce an analytically solvable model for SHG 
in localized waves, we follow Ref. ITS and consider a fun- 
damental frequency (FF) wave propagating along the Z- 
direction in a linear slab waveguide, as shown in Fig. [I]. 
We assume that the interfaces (or defect layers) possess 
a quadratic nonlinear response, so that a FF wave can 
parametrically couple to its second harmonic (SH) wave. 
The coupled-mode equations for the complex envelope 
functions Ej(x, Z) (we use j=l,2 for FF and SH, respec- 
tively) can be written in the form, 



dE d 2 E 

+ £>i^r + si(x)Ex + T 1 {x)ElE 2 = 0, 



dZ 
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l ~dz +D2 - 
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+ e 2 (x)E 2 +T 2 (x)E 2 = 0, 



du d 2 u 
l J- Z + dx- 2+Vl[x)u 



+ S (x — x n ) {(3\u + u*v) = 0, 



.dv d 2 v 
l Yz +a dx 27 + " 2{x)v 



(2) 



+ J2 S ( X - X n) (02V + U 2 ) = 0. 



At this point, it is important to note that our equa- 
tions [Eqs. (Q) or (g)] describe the beam evolution in 
the framework of the so-called parabolic approximation, 
valid for the rays propagating mainly along the Z di- 
rection. In other words, the characteristic length of the 
beam distortion due to diffraction and refraction along 
the Z axis should be much larger than the beam width 
in the transverse direction x. This leads to the require- 
ment of a shallow grating, \eoj(x) — eoj\ <SC £oj- On the 
other hand, to make the parametric interaction effective, 
the mismatch between the phase velocities of the FF and 
SH waves should be small, 2eoi — £02 • Then, the ratio of 
the diffraction coefficients is approximately a — 1/2, and 
we use this value in the numerical simulations presented 
below. 

I |] linear medium 
I \ — interface 




FIG. 1. Schematic diagram of an array of thin \^ layers 
embedded into a linear slab waveguide, with indication of the 
direction of the input laser beam relative to the interfaces. 



where Dj are the diffraction coefficients (Dj > 0). In 
the approximation of infinitely thin interface layers (valid 
when the width of each layer is much smaller than the 
characteristic transverse scale of the FF and SH wave en- 
velopes), we take £j(x) — £oj{x) + J2 n K j$( x ~ x n) an d 
Tj(x) = ^2 n ^j5{x — x n ), where x n denotes the position 
of the n-th nonlinear interface, jj are the nonlinearity 
coefficients, £oj(x) and Kj account for the phase velocity 
differences in bulk and interface materials. 

In order to reduce the number of physical param- 
eters, we normalize Eqs. (|l|) as follows: E\(Z) = 
u(z)e l ^ z /^Wn, E 2 (Z) = v{z)e 2 y^ z hx, a = 
D 2 /Di, Vj(x) = (e j{x) - j £ i)/-Di, % is the average 
value of eoj(x), and (3j = Kj/D\, where z = Z/D\ is mea- 
sured in units of D\ . Then, the dimensionless equations 
become 



For spatially localized solutions, the system (|2|) con- 
serves the Hamiltonian 



H 





du 


1 


a 


dv 


r{ 

1 —00 1 


dx 


H 


' 2 


dx 



vi{x)\u\ —\v\ 



■y~]6(x- X n ) 



/3i|u| 2 + ^|^| 2 + Re( U V) 



dx 



and the total power P = f_°° (\u\ 2 + |i>| 2 ) dx. 



III. PERIODIC LAYERED STRUCTURES 
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A. General formalism 



To develop a general formalism for describing station- 
ary, spatially localized modes, we consider an infinite sys- 
tem of uniformly spaced nonlinear interfaces located pe- 
riodically at the positions x n = nh, which separate iden- 
tical linear layers, i/j (x + nh) = i/j (x). Example of such 
a structure is shown in Fig. ^. This type of a one- 
dimensional (ID) nonlinear photonic crystal (NPC) can 
be used to prohibit light propagation along the transverse 
x axis under certain conditions, resulting in field local- 
ization, and it resembles the operation of the so-called 
photonic crystal fiber in the linear regime |l4f| . The fun- 
damental properties of NPC can be understood by study- 
ing nonlinear localized modes. Such modes appear in the 
frequency gap of the layered linear structure, and they 
can be called discrete gap solitons or, in a broader con- 
text, intrinsic localized modes [01. 




FIG. 2. The slab waveguide from Fig. [j] viewed from the 
front with indication of the substructure and notation of the 
identical linear media (photonic crystal with M sub-layers 
with different refractive indices) in between the x' 2 ' inter- 
faces. The separation between the nonlinear interfaces is h. 
A localized mode centered at x~0 is also illustrated. 



We search for solutions of Eqs. (g) in the form: 

u (x, z) = d (x) e iXlZ , v (x, z) = c 2 (x) e lX2Z , 

where Cj (x) describe the transverse profiles of FF and 
SH waves, respectively, and the real propagation con- 
stants Xj satisfy the phase-matching condition A2 = 2Ai . 
Then, as the FF and SH fields do not interact in a linear 
bulk medium, each of them can be presented as a com- 
position of two linear eigenmodes, corresponding to pairs 
of counter-propagating waves with opposite wavevector 
components along the x axis. Regardless the internal 
structure of a linear layer between nonlinear interfaces, 
the wave amplitudes at its boundaries can be related by 
a transfer matrix: 



,(«+!-) 



n+1,-) 



>>+) 



,+) 



(3) 



Here indices ± stand for the waves on the right and on 



the left of a given nonlinear layer, respectively, and the 
amplitudes of counter-propagating waves are denoted by 
Oj and bj, as shown in Fig. ||. 

In order to calculate the dependence of the matrix ele- 
ments on the propagation constants, we should solve the 
corresponding linear problem. For practical applications, 
NPC can be produced by embedding nonlinear layers in 
an otherwise linear Bragg grating structure, see, e.g., 
Ref. [^2[ . Thus, we assume that each linear layer con- 
sists of several sub-layers with a constant refractive in- 



dex, see Fig. g, i.e. Vj(x) 



for 



< x < x 



where m is used to number the sub-layer inside a lin- 
ear slice (1 < m < M), and we have, by definition, 
x n .i = ^n-i. M+i = x n = nh. Then, the field in each 
of the sub-layers can be written as: 

Cj (x) = a ^ e -HM— n ^) + ^"V^ (4) 



By definition, the amplitudes from Eq. (H) are a 



(n,l) (n+1,-) 



ing for b 



0,±) 



(n,M+l) 
2j , 
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with similar relations hold- 



a 



(n,+) 



b 



(n,+) 



a 



(n+1,-) 



b 



(n+1,-) 



(.1) 



(n+1) 



FIG. 3. Transfer matrix relating the amplitudes at the 
boundaries of a linear layer, independently for the FF and 
SH components, see Eq. (^). 



Transverse wavenumbers can be calculated using the 
dispersion relations, Ai = /xf m + v\^ n and A2 = cr/x| m + 
i/2, m , together with the phase matching condition A2 = 
2Ai. From these expressions we conclude that, in general, 
for given Xj, the waves can be either localized (i.e., fij. m 
is real) or propagating (i.e., /Uj, m is imaginary). Then, we 
notice that the waves in a multi-layered linear medium 
can be determined as linear eigenmodes localized at the 
sub-layer boundaries. First, from Eq. (Q) it follows that 
the variation of the field amplitude through the m-th sub- 
layer of the width h m = x n , m +i — x n m is characterized 
by the following transfer matrix: 



p 



0. 







(5) 



Second, the variation of the wave amplitudes at the 
boundary can be calculated by applying the field continu- 
ity conditions following from the model (|^). Specifically, 
we equate the amplitudes Cj and their derivatives dcj/dx 
on both sides of the interface, and find the transfer ma- 
trix accounting for the field localization at a boundary 
separating the m-th and (m + l)-th sub-layers, 
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(6) 



where P t ' m = Pj. m / Pj.m+i- To simplify the further analy- 
sis, we include all linear properties in the transfer matrix. 
Then, Pj y M+i = due to periodicity of the underlying 
grating, and the linear (^-response is characterized by the 
following matrix: 



We note that some particular cases of the system (g) 
have been earlier discussed (but, in fact, never derived 
in a consistent manner) in the analysis of the nonlin- 
ear interface dynamics under the condition of Fermi res- 
onance ]l5| , |l6lh arrays of weakly interacting quadratic 
waveguides |17| pof , and beam propagation in nonlinear 
lattices il. 



T, 



where ty = 



1 



l 8 



(7) 



and m is the index of the linear 



layer with the delta-interface (in our case, M + 1). The 
total transfer matrix can be then found as a product, 



T UM+l) T (j,M) T (j,M) 



T(j,i)T(j,i) 

' r p 



(8) 



After calculating the matrix elements, we express the 
field in terms of the mode amplitudes defined at the 
nonlinear interfaces, = Cj (x n ), combining the re- 
lation (H) with the continuity conditions at the layers, 



» 



(n,M+l) 



,(n,M+l) t i , 

b\ . Integrat- 



7n,l) , (n,l) 

°j - a j + °3 -3 1 "3 

ing Eqs. (g) over small intervals including the nonlinear 
interfaces at the positions x n (note that we exclude lin- 
ear responses as they are already accounted for in the 
transfer matrices), we derive the so-called discrete x^ 
equations: 



VlUn + (U n -1 + U n+1 ) + XlU*V n = 0, 

mv n + (K-i + v n+1 ) + = 0, 



(9) 



where U n = -\/l^i^2] and V n — c£ are the nor- 
malized FF and SH amplitudes at the n-th nonlinear 
layer. Parameters Xj and rjj are defined by the matrix 
elements, 



% = -Tr(T«), x, =sign(0), 



(10) 



where 



T U) = T 



(i) 



• T, 



(3) 



[2,1] 



T, 



T, 



(3) 

[2,2] 



0) 
1,2] 



T, 



(3) 
2,2]' 



(11) 



B. Discrete gap solitons 

Now we use Eqs. (^) to find the stationary localized 
modes of nonlinear superlattices, or discrete gap solitons. 
A similar problem was earlier analyzed in Ref . [ fi"5| in the 
so-called continuum limit, when the modes become wide 
and are effectively supported by many interfaces, with 
the excitation profiles approaching those of quadratic 
solitons Q . On the other hand, it has been demonstrated 
that discrete states in a closed system of few interfaces 
can have different topologies and possess quite peculiar 
properties jL(|[l7). For the case at hand, when the num- 
ber of interfaces is infinite (e.g. much larger than the 
characteristic mode width), the different types of highly 
localized waves have been identified pjjfl , and their pro- 
files were described by approximate analytical solutions. 
However, until now the transitional case of moderately 
localized modes has not been addressed. Thus, we de- 
velop a more complete analytical description of discrete 
gap solitons, which can predict their properties in all the 
parameter regions. 

In order to find approximate solutions for highly lo- 
calized modes, we use the variational method. First, we 
have to choose the trial functions. We use the fact that in 
the high localization limit the tails are almost linear, so 
that the amplitudes decay nearly exponentially. Then, 
we introduce two sets of trial functions to account for 
different topologies |l8|]: odd modes, when the center of 
symmetry is located at a layer, 



(12) 



and even modes, when the center of symmetry is located 
between two neighboring layers, 



(o) _ J C/ 4 ,l| e- plH , n> 0, 



U, 



V, 



U ts l " + 1| e - pi l" +1 l, n < 0, 
Vb4"'e" P21 ™ 1 , n> 0, 

-p 2 |n+l| ) n < q 



(13) 



VnS 



0*2 



|n+l| 



Hereafter we use the notation T, , to denote the matrix 

[n.m\ 

element in the row n and column m of the matrix T^K 
It is easy to verify that the parameters rjj and normal- 
ization coefficients £j are real for any (real) propagation 
constants Xj. A proof of this fact, together with discus- 
sion of some other properties of the constructed transfer 
matrices, is given in Appendix A. 



Here the parameters Sj — ±1 are introduced to describe 
unstaggered and staggered profiles, and t = ±1 to pro- 
duce either untwisted or twisted modes (for the signs "+" 
or "— ", respectively). 

After selecting the mode topology and fixing the values 
of Sj and t, the unknown values Uq, Vq, Pj (pj > 0) are 
determined by minimizing the Lagrangian corresponding 
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to Eqs. (|9|). Details of these calculations will be presented 
elsewhere p2| , here we give a brief summary of the main 
results, and discuss their physical consequences. 

As the system @ possesses the symmetries (i) \i — ¥ 
~Xi, s i -* si, Vi ~> -Vi, and (ii) \2 -> ~Xa and 
— * — V n , we consider, without a lack of generality, the 
case Xj — 1- The analysis shows that localized solutions 
exist only if s 1.771 < —2. The latter condition means that 
the FF component is unstaggered for 771 < —2, and stag- 
gered, otherwise (771 > 2). Similarly, we consider only 
the case 1 772 1 > 2, as for other values the localized solu- 
tions are unstable due to resonant interaction with linear 
waves Q. 



be found as V n ~ L^/772, resulting in an effectively cubic 
nonlinearity for the FF wave (lj] (see also [Q). In the 
intermediate case, a transition of the SH profile between 
staggered and unstaggered topologies should be observed. 
Indeed, our variational calculations predict that the SH 
is staggered for 2 < 772 < 7722, and unstaggered for 
772 > 7722- Here, the critical parameter value depends 
on 771: for odd modes, it is found from the quadratic 
equation, ^7722 + 1 / ^7722 = |?7i |, and a similar relation 
holds for even modes V^22 + 1 + l/Vte + 1 = The 
topology does not change sharply as the parameter 772 
crosses the critical value 7722 (771), but numerical calcula- 
tions confirm that such a transition occurs in a region 
close to the separation line. 
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FIG. 4. Different types of odd two-frequency localized 
modes (numerical solutions: FF - triangles, SH - squares; 
variational profiles: FF - dashed line, SH - solid), (a) stag- 
gered FF and unstaggered SH (771 = 2.4, 772 = 4.5), (b) both 
components staggered (771 = 2.4, r/2 =3), (c) both unstag- 
gered (771 = —2.4, 772 = —3.5). 

Analyzing the linear problem, it is straightforward to 
see that the SH mode can be staggered only if 772 > 2. 
Then, we distinguish between two limits: (i) a strongly 
localized FF mode (771 ^> 2), when the SH consists of 
staggered linear tails [ pj| |, and (ii) the cascading limit 
(772 ^S> 2), when the SH profile is unstaggered and can 
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FIG. 5. Characteristic examples of even two- frequency lo- 
calized modes with untwisted FF component; the notations 
are the same as in Fig. ^. Parameter values are (a) 771 = 2.45, 
772 = 4.5, (b) 771 = 2.45, 7/2 = 2.5, (c) 771 = -2.45, 772 = -3. 

We performed numerical analysis and found that the 
approximate variational solutions provide close match- 
ing for the highly localized profiles, i.e. for relatively 
large values of Examples of odd and even modes 

are presented in Figs. ^ and ||, respectively. We see that 
the profiles of staggered modes, supported by only a few 
interfaces, are described very accurately. However, for 
wider modes the deviations between the exact numeri- 
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cal and variational solutions are more pronounced, see 
Figs. |](c) and ||(c). The limitation of the variational 
solution is due to the specific choice of the trial func- 
tions (fL2|) , (|l3|) , which arc not suitable for description of 
"moderately" localized waves with smoother profiles. 

It can be demonstrated that in the continuum limit 
(jjj — > — 2") untwisted modes acquire the profiles of 
quadratic solitons , that can be well approximated with 
the sech-type functions [p3[| , A special quasi- continuous 
approach, which allows to determine the mode profiles as 
a soliton bound state can also be developed. The result- 
ing approximate solutions provide very good estimates, 
and they are very useful for understanding the mode scal- 
ing properties, i.e. a change from broad solitons to nar- 
row highly localized states. On the other hand, it should 
be noted that the twisted modes do not exist close to 
the continuum limit, because their profiles are intrinsi- 
cally discrete due to a sharp amplitude change between 
the layers at the mode center. A comprehensive descrip- 
tion of these results goes beyond the scope of the present 
paper, and will be presented elsewhere |22]. 



IV. A SINGLE x (2) NONLINEAR LAYER 
EMBEDDED IN A PERIODIC STRUCTURE 



(14) 



Here JV^) is the transfer matrix of one linear seg- 
ment in a periodic lattice, starting on the right (+) or 
the left (— ) side of a nonlinear interface, and it can 
be calculated using Eqs. (Bj)-(^). Then, we notice that 
the determinants are: \\To ™' ^11 = 1, 



IT, 



0',m,±)i 



J^ m+1 , and thus ||T^±)|| = 1. Using 
this relation, we solve Eq. (jl^) and find the values of the 
reflection coefficients: 



_± _ T (j,±) 



rU,±) 



[2,2] 



(15) 



where 



(vt/2) 



1-Wl- 2/77 



It can be proved that rfj are real (see Appendix A). Then, 

we see that the modes decay at infinity, i.e. Ir^'^l < 1, 
only if \f]j\ > 2. This condition, which has already been 
mentioned in Sec. HI, defines a band gap, when the wave 
propagation is prohibited for a specific range of the prop- 
agation constant Aj (see an example in Fig. ^|). 

Band structure \ 



Let us now consider a special situation, when there 
exists only a single nonlinear layer located at xq = em- 
bedded into a linear grating. A similar problem has been 
considered in our recent paper for the case when 
the linear media on both sides of the nonlinear inter- 
face are uniform. Here, we generalize those results to 
the case of non-uniform linear media, considering a non- 
linear layer with quadratic nonlinearity embedded into a 
periodic structure. Therefore, we modify the model (|^) 
as follows: (i) nonlinear coupling terms are taken into 
account only for n = 0, and (ii) a linear response of the 
central layer is assumed to be different from that of other 
layers, i.e. we change (3j — ► aj + j3j at n — 0. In connec- 
tion with the previous problem (see Sees. || and III), the 
case with aj = corresponds to the limit of a highly lo- 
calized mode, when the mode width is much smaller than 
the distance in between the nonlinear layers (see Fig. |^). 

Following the general approach outlined in the previ- 
ous section, we first analyze the linear properties. We 
also use the similar notations, but for a single nonlin- 
ear interface we omit the index n = 0. From the theory 
of linear structures, it follows that a link between the 
linear-wave amplitudes can be characterized by the re- 
flection coefficients, r+ = b^/a^ and rj — aj /bj . Wc 
do not assume that the linear structure is symmetric [in 
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FIG. 6. The band structure for (a) FF and (b) SH modes, 
and (c) combined band gap (additional narrow gaps in FF and 
SH spectra exist for smaller Ai, not shown). Areas without 
shading correspond to band gaps. Parameters are h± — 0.6, 



ui,+ = 0, fi. 



-6, V2, 



-8, v 2 , 



-13, and fl a 



0. 



By applying the continuity conditions, Cj = a - +b ■ , we 
obtain the expressions for the amplitudes at the nonlinear 



layer: \c\ | 2 = a\ai, C2 = — a\, where a\ — a\ — — ( 



Cl2 



Ct2 



are the effective transverse wave numbers at the right 
(+) and left (— ) boundaries of the nonlinear layer; their 
values are determined by the dispersion relation of the pe- 
riodic linear gratings. As has been demonstrated in the 
general Vj(x) ^ Vj(-x)] and denote with + and - the previous section, the wave numbers ^^_the linear trans- 
wave characteristics at the right and left boundaries of 
the nonlinear layer. If the linear structure is periodic, 
the coefficients can be found by solving the following 
eigenvalue problems: 



51) , the reflection 



fer matrices, and, according to Eq 

which, in turn, are related by the phase-matching condi- 
tion A2 = 2Ai. Thus, for fixed physical characteristics, 



coefficients rj~ depend on the propagation constants A^-, 
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the localized modes constitute a one-parameter family, 
and we choose Ai as a free parameter. We note that the 
coefficients £ x are real in the band gap (see Appendix B). 




j i i i i i i i i i i i i i 

-4 4 8 12 16 20 




FIG. 7. Three types of the power dependences for 
two-color stationary localized modes, for different linear mis- 
matches: (a) Qi = — 1 and a?2 = — 1; (b) a?i = 10 and ai = 8; 
(c) «i = 4 and a?2 = 6. Solid - stable, dashed - unstable, 
and dotted - oscillatory unstable. Parameters of the linear 
structure correspond to those in Fig. 

In order to illustrate the features of nonlinear modes, 
we consider a structure similar to that used in exper- 
iments on the SH generation in localized modes [ jl3| . 
In the experimental setup, a periodic linear grating was 
built of two materials with different refractive indices, 
characterized by the parameters Vj ±, with correspond- 
ing finite widths of the layers, h±. The nonlinear inter- 
face was created by cutting the grating in two parts and 
coating the interface to enhance the effective quadratic 
nonlinearity. Characteristics of the defect layer were con- 
trolled by adjusting the gap. An example of the band gap 



for such a photonic structure is presented in Fig. |6|. Note 
that because the linear structures on either sides of the 
nonlinear layer are chosen to be the same [up to a con- 
stant shift, in our notation i/j i —(x) — Vj t +(x + h+)], the 
corresponding gaps coincide. 

Many of the properties of the localized states can 
be understood by analyzing the power diagram P(Ai). 
Characteristic examples of such a dependence are shown 
in Figs. f?|(a,b,c). Similar to the case of homogeneous 
linear media p3|, there always exists a branch in a pa- 
rameter region unbounded from above for Ai larger than 
some critical value. Quite remarkably, the corresponding 
mode properties are very similar to those of quadratic 
solitons j7j. In particular, for large values of the propa- 
gation constant Ai the power dependence P (Ai) always 
has a positive slope, while for smaller Ai the slope can be 
negative, resulting in bistability. Such a case is demon- 
strated in Fig. @(a). We found that stability of the cor- 
responding modes can be determined by the Vakhitov- 
Kolokolov criterion, i.e. the localized modes are stable 
provided dP/dXi > 0, and unstable, otherwise. 




15 30 



z 

FIG. 8. (a) Generation of a stable localized mode from 
an input FF beam with a Gaussian transverse profile (ini- 
tial power P ~ 8); (b) Switching from a perturbed unstable 
(corresponding to Ai = —I) to a stable mode. Left: the field 
intensities of the FF (dashed) and SH (solid) components at 
the interface. Right: evolution of the SH field. Parameters 
correspond to those in Fig. ^(c). 

On the other hand, the spectrum of a linear periodic 
structure consists of several bands. Moreover, even in- 
side a band gap the modes cannot exist if the condition 
ct\oi2 > is not satisfied. Thus, in a sharp contrast 
with two-color parametric solitons in homogeneous me- 
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dia, other branches can appear for smaller values of Ai, 
and we found that the localized mode properties can be 
very different compared to the modes corresponding to 
the right branch. For example, for Figs. 0(b,c), the left 
branches at higher intensities correspond to smaller Ai 
and wider profiles. This happens because the SH am- 
plitude is negative, C2 < 0, which results in effectively 
self-defocusing nonlinear response. It is interesting to 
note that similar types of power dependencies occur for 
a layer possessing a self-defocusing Kerr-type nonlinear- 
ity |24j , and the corresponding localized modes are stable. 
We have performed a linear stability analysis and have 
found that, similar to the case of the Kerr-type nonlin- 
earity, the Vakhitov-Kolokolov type criterion cannot be 
applied to determine stability of such modes existing in 
an effectively defocusing medium. However, our calcu- 
lations demonstrate that parametric resonances can lead 
to oscillatory instability, as illustrated in Figs. |?](b,c). A 
comprehensive discussion of the mode stability will be 
presented elsewhere. 

The two-color modes can be generated by launching a 
localized FF wave at the interface, as shown in Fig. ||(a). 
We have also studied the evolution of perturbed unstable 
modes, which can evolve toward a stable state. Example 
in Fig. ||(b) demonstrates the development of an oscil- 
latory instability with subsequent switching to a stable 
branch. 



We believe our results are important, on one hand, for 
the theory of nonlinear photonic crystals where nonlin- 
earities appear due to phase-matched harmonic genera- 
tion, and, on the other hand, for creating tunable band 
gap materials where gaps could be opened or closed de- 
pending on the input intensity. 
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APPENDIX A: TRANSFER MATRIX 
PROPERTIES 



V. CONCLUSIONS 



We have developed a general formalism for analyzing 
spatially localized nonlinear modes (discrete gap solitons) 
in periodic photonic structures with embedded quadratic 
(or x^ 2 -*) nonlinear interfaces — nonlinear quadratic su- 
perlattices. Our approach can be applied to different 
types of periodic linear media with isolated or peri- 
odic nonlinear interfaces, where nonlinearity can sup- 
port parametric interaction and generation of the second- 
harmonic field. In the case of a nonlinear superlattice, 
i.e. periodically spaced layers with quadratic nonlinear 
response, we have derived an effective discrete model. 
For such a periodic structure, we have found discrete gap 
solitons of different topologies in the form of fundamen- 
tal and second-harmonic fields coupled parametrically at 
the nonlinear interfaces. 

For a single nonlinear layer embedded in a linear peri- 
odic medium, we have described a novel type of nonlin- 
ear localized defect mode — a two-color photonic mode. 
Some of the properties of these two-color localized modes, 
such as stability, generation, and switching, have been 
shown to be remarkably similar with those of quadratic 
parametric solitons in homogeneous media. However, we 
have revealed a number of specific properties of such 
modes, and also demonstrated their generation from a 
localized beam of the fundamental frequency, as well as 
switching from an unstable to a stable state. 



In order to demonstrate some features of linear modes 
existing in periodic structures, we consider the proper- 
ties of the corresponding transfer matrices. First, we see 
that a matrix describing changes of the wave amplitudes 
at the boundary between linear layers, T r , depends only 
on the wave numbers on either side of the interface, as 
follows from Eq. (||). It is also easy to check that the 
following relation holds: 

T r ^' ^ = -^r(M7,m+l s M?,m) 

T r {^j,m+1, ft) T r (ji, flj^ m ), 

where ft is arbitrary, and we assume that it is 
real. Let us introduce a new matrix Tw' = 
T T (ft, fij t M+i)T^T T (fij t i, ft), and use Eq. (|l^) to present 
it in a special form: 

Here f^ n) = T T (ft, Mj , „)T^ n) T r (^- n , ft) and Tp' n) = 

T T (ft, /j,j jn )Tg n 'T T ((j,j tn , ft). It can be verified by direct 
substitution that these matrices are real (we use the fact 
that, for stationary modes, /j,j iTl are real or purely imagi- 
nary), and therefore T^) is real as well. Finally, we find 

that Tr(T^) = Tr(T and T s (i) /Mj,i = T^/ft, which 
proves that coefficients l~j and rjj in Eq. (^) are real. 
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APPENDIX B: REFLECTION COEFFICIENTS 
FOR MODES IN A BAND GAP 



Although it is possible to extend the technique pre- 
sented in Appendix A to prove the properties of the re- 
flection coefficients in the case of an infinite linear grating 
(introduced in Sec. IV), here we employ a different ap- 
proach. We note that for stationary waves in a band gap 
there should be no energy flow along the x axis. The 
corresponding restrictions follows from Eq. (H): 



Im 



du* 
dx 



0, Im 



dv* 

) 

dx 



0. 



(17) 



These conditions are satisfied when the amplitudes of the 



counter-propagating waves coincide, i.e. 



1 if \x 



± 

3 

is imaginary. On the other hand, for the layers with real 
//^ the linear modes should be in phase, i.e. Im(r^) = 0. 
Then, it immediately follows that in a band gap the co- 
efficients Cf are real. 
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